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The maintenance of the global H 2 SO 4 clouds on Venus requires volcanism to re- 
plenish the atmospheric SO 2 which is continually being removed from the atmosphere 
by reaction with calcium minerals on the surface of Venus . 1-4 The first laboratory 
measurements of the rate of one such reaction, between SO2 and calcite ( CaCO 3 ) to 
form anhydrite ( CaSO 4 ), are reported here. If the rate of this reaction is represen- 
tative of the SO 2 reaction rate at the Venus surface, then we estimate that all S0 2 
in the Venus atmosphere (and thus the foSO* clouds) will be removed in 1.9 million 
years unless the lost SO 2 is replenished by volcanism. The required rate of volcanism 
ranges from about 0.4 to about 11 km 3 of magma erupted per year, depending on the 
assumed sulfur content of the erupted material. If this material has the same compo- 
sition as the Venus surface at the Venera 13, 14 and Vega 2 landing sites 5,6 then the 
required rate of volcanism is about 1 km 3 per year. This independent geochemically 
estimated rate can be used to determine if either (or neither) of the two discordant 
(2 km 3 yr -1 vs. 200-300 km 3 yr -1 ) geophysically estimated rates 7-9 is correct. The 
geochemically estimated rate also suggests that Venus is less volcanically active than 
the Earth. 

The experiments utilized calcite, which is a net sink for SO 2 by incorporating it into the Venus 
crust via the net reaction 1-4 

CaCOz (calcite) 4- S0 2 — ♦ CaSC >4 (anhydrite) +CO. (l) 

Calcite is also predicted by several authors to be present on the Venus surface because it ap- 
pears capable of buffering the C0 2 abundance on Venus via the Urey reaction (CaCO^ + 5i0 2 = 
CaSiOs + C0 2 ). 10-13 These predictions concerning possible roles for calcite on Venus, the ex- 
emplary nature of reaction (1), and the availability of high purity CaC0 3 in several forms led us 
to study reaction (1) as a paradigm for S0 2 removal. However, as Figure 1 illustrates, S0 2 also 
reacts with other calcium minerals (e.g., diopside CaMgS^Os), which are believed to be present 
on the surface of Venus. 14,15 The kinetics for these alternative reactions are also being measured 
and will be reported later. 

Reaction (1) was studied by heating clear calcite crystals (Iceland spar) of known impurity 
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content (~ 350 ppm total impurities by inductively coupled plasma emission spectroscopy), weight, 
and surface area in S0 2 -bearing gas streams for varying time periods at three different tempera- 
tures. Experiments were done at temperatures of 600-850° C and ambient atmospheric pressures 
with analyzed S0 2 - C0 2 gas mixtures nominally containing 1% (by volume) S0 2 . Thus, the 
SO 2 molecular number density is about the same as at the surface of Venus where the atmospheric 
pressure is about 100 times larger but the S0 2 volume mixing ratio is about 100 times smaller. 16 
The mixtures also typically contained < 1 ppm 0 2 and < 10 ppm if 2 0 as impurities. Standard 
techniques were used to monitor and control gas flow and temperature (to better than ± 5°C). 
The rate of reaction (1) was determined by three independent methods: (a) measuring the weight 
gain at the end of reaction, (b) using scanning electron microscopy (SEM) to measure the porosity 
and thickness of the reacted surface layers from micrographs of external and fracture surfaces (e.g., 
Fig. 1), and (c) analyses of reacted samples by ion chromatography. SEM examination of 

many samples (in combination with X-ray line scans and energy dispersive spectroscopy (EDS) on 
the SEM) showed that CaSO A layers covered all external surfaces of the reacted samples. X-ray 
powder patterns of samples and standards confirmed that the layers are anhydrite (in particular, 
there is no evidence for calcium sulfite) and that unreacted calcite is the only other phase present 
in the samples. The ion chromatography also verified the absence of anion impurities at levels > 
10 ppm. 

The experimental run conditions and rate data are summarized in Table 1 and Figure 2. Rate 
data using the weight gain, SEM and SO\~ analysis methods for the same sample agree within the 
combined experimental uncertainties (la) shown in Table 1. A weighted linear least squares fit 
to the data gave the equation for the rate R = io 19 - 64 l ±0 - 28 ) exp (-15,248 (± 2970)/T) molecules 
cm -2 s _1 . The corresponding activation energy is 126.8 ± 5.4 kJ/mole. The rate of reaction (1) 
was extrapolated downward to Venus surface temperatures (~ 660-750 K) using this equation and 
the Pioneer Venus radar altimetry data 17 and atmospheric (P, T) profile 18 to take the altitude 
dependence of the rate into account. The derived global mean S0 2 reaction rate is ~ 4.6 x 10 10 
molecules cm -2 s _1 . 

Chemical analyses of the Venus surface by the Venera 13, 14 and Vega 2 spacecraft 5,6 give 
a CaO content (by mass) of 7.90% (weighted mean). Assuming the experimental rate is repre- 
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sentative of the SO 2 depletion rate via reaction with calcium minerals on the Venus surface, the 
observed SO 2 column density (2.2 x 10 23 molecules cm -2 ) 16 would be removed from the Venus 
atmosphere in ~ 1.9 x 10® years in the absence of a comparable sulfur source. However, main- 
tenance of the global H 2 SO 4 clouds, which are formed by the UV sunlight powered conversion 
of SO 2 into H 2 SO 4 cloud particles, 19 requires a comparable sulfur source. This source must be 
endogenic because infalling cosmic material cannot supply enough sulfur. The calculated SO 2 
removal rate on Venus is equivalent to ~ 2.8 x 10 13 grams sulfur per year. However, the measured 
terrestrial flux of infalling material, which is roughly applicable to Venus because of similar sizes 
and masses, is only — 7.8 x 10 10 grams per year, 20 or only ~ 4.6 x 10 9 grams sulfur per year 
assuming the most sulfur-rich Cl chondrite composition. 21 

The most plausible endogenic source is volcanism, which has occurred on Venus in the past, 22 
and which may have led to increased SO 2 levels above the Venus cloud tops observed by the PV 
orbiter. 23,24 The rate of volcanism required to balance SO 2 depletion via reactions with calcium 
minerals on the Venus surface depends on the sulfur content of the erupted material (gas and 
magma). If this material has an overall S/Si mass ratio of 0.03, the weighted mean of the Venera 
13, 14, and Vega 2 surface analyses, 6,6 the corresponding rate of volcanism is ~ 1 km 3 yr' 1 . Figure 
3 illustrates that this rate is about 20 times smaller than the terrestrial value. Two other plausible 
models for the S/Si ratio of the erupted material on Venus are also shown: (a) S/Si ~ 0.1, as in 
ordinary chondrites, 21 giving a rate of ~ 0.4 km 3 yr -1 , and (b) S/Si ~ 0.004 as in the Earth’s 
crust, 25 , giving a rate of ~ 11 km 3 yr -1 . All rates will scale as e -1 , where e is the S degassing 
efficiency for erupted material (O < e < 1). However e is likely to be close to unity because S gases 
are expelled directly into the atmosphere and hot, reactive S-bearing magmas will rapidly lose S 
to the atmosphere. Our experiments on a possible magma S-bearing mineral pyrite FeS 2 being 
chemically weathered by CO 2 indicate complete S loss in hours at sub-magmatic temperatures (~ 
980° C). 

Figure 3 also illustrates that widely discordant volcanism rates have been estimated by two 
different geophysical methods: (a) rates of ~ 200-300 km 3 yr -1 by scaling terrestrial heat pro- 
duction to Venus and transporting the outward heat flux solely by volcanism 8,9 and (b) an upper 
limit of 2 km 3 yr -1 by analysis of impact crater areal densities on Venera 15, 16 radar images 7 . 
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The present geochemically derived volcanism rate is independent of these two methods and in 
principal can be used to test if either (or neither) is correct. In particular, the values computed by 
geophysical method (a) do not agree with our geochemical values. Finally, this work also suggests 
that on average Venus is less volcanically active than the Earth. 

We thank W. Zamechek for chemical analyses, and L. Rigione, S. Rowell, and A. Petropoulos 
for assistance in the laboratory. This work was supported by the NASA Innovative Research 
Program. The MIT UROP program also partially supported S.R. and A.P. 
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TABLE 1 


Rate Data Used to Calculate the Activation Energy 


Samole 

T(°C) a 

Duration 

(Hours! 

Rate xlO -12 
(molecules cm' 2 s - 

R44-1 

602 

95 

1.67 

± 0.19° 

R45-1 

602 

314 

0.85 

± 0.77 




1.05 

± 0.12 c 

R45-3 

602 

314 

1.01 

± 0.11 c 

R45-4 

602 

314 

1.20 

± 0.51 




0.68 

± 0.48 d 

R36-1 

746 

96 

8.43 

± 2.93 




13.6 

± 6.0 d 

R36-2 

746 

96 

23.2 

± 6.2 




19.9 

± 11. O'* 

R50-1 

747 

192 

9.19 

± 1.60 

R50-3 

747 

192 

10.5 

± 7.8 

R51-1 

748 

96 

13.2 

± 4.0 

R51-2 

748 

96 

10.3 

± 5.6 

R43-1 

848 

91 

60.4 

± 6.8 




57.1 

± 31.7 a 

R43-2 

848 

91 

53.2 

± 7.0 




67.2 

± 43. 4 d 

R43-3 

848 

91 

58.1 

± 7.0 




48.4 

± 12. 1^ 


Notes to Table 1 

a The estimated temperature uncertainty is ± 5°C. 

6 The number of SO 2 molecules reacted per cm 2 surface area per second. The rate is from 
weight gain measurements unless noted otherwise. The 1 a uncertainties are also listed. 

c The rate from ion chromatography analyses of SO 2- content of reacted sample. 

d The rate from scanning electron microscope (SEM) measurements of anhydrite layer thickness 
and porosity. 



Figure Captions 


Figure 1. Scanning electron micrographs of chemically reacted samples, a, Fracture surface showing 
a CaSOi (anhydrite) layer (A) completely covering a partially reacted calcite crystal (C) (850° C 
for 192 hours) b, Top surface showing CaSO 4 (A) growing on a partially reacted diopside crystal 
(D) (833° C for 48 hours). 

Figure 2. Arrhenius plot of rate data from Table 1 vs. 1/T. The line is a weighted linear least 
squares fit. Rates determined from independent methods (wt. gain, SEM, SO^~ analyses) for the 
same sample agree within the 1 a experimental uncertainties. 

Figure 3. Geochemical and geophysical estimates of the rate of volcanism on Venus are compared 
to the terrestrial rate. Three plausible models for the S/Si ratio of erupted material on Venus 
yield rates of ~ 0.4 to ~ 11 km 3 yr. -1 . Adopting the Venus surface S/Si ratio from Venera 13, 
14, and Vega 2 analyses 5,6 yields a rate of ~ 1 km 3 yr. -1 , or about 20 times less than the rate of 
volcanism on Earth. 
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